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Lemaitre-Robertson- Walker universe admits one, and only one, gauge-invariant variable which de- 
scribes the perturbation to the energy density and which becomes equal to the usual Newtonian 
energy density in the non-relativistic limit. The same holds true for the perturbation to the particle 
number density. Using these two new variables, a manifestly covariant and gauge-invariant cosmo- 
logical perturbation theory, adapted to non-barotropic equations of state for the pressure, has been 
developed. The new perturbation theory is free from metric gradients. 

Perturbations in the total energy density are gravitationally coupled to perturbations in the 
particle number density, irrespective of the nature of the particles. Small-scale perturbations in the 
radiation-dominated era oscillate with an increasing amplitude. 

After decoupling of matter and radiation density perturbations evolve diabatically, i.e., they 
exchange heat with their environment. This heat loss of a perturbation may enhance the growth 
rate of its mass sufficiently to explain stellar formation in the early universe, a phenomenon not 
understood, as yet, without the additional assumption of the existence of Cold Dark Matter. This 
theoretical observation is the main result of this article. 
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I. INTRODUCTION 

Since measurements of the fundamental parameters of our universe are very precise today, cosmology is nowadays 
a mature branch of astrophysics. Despite advances in observational as well as theoretical cosmology, there is, as 
yet, no consensus about which cosmological perturbation theory describes unequivocally the evolution of density 
perturbations. 

A new approach to cosmological density perturbations is presented and evolution equations are given which are 
adapted to the situation that the equations of state are more realistic. The issue of metric gradients [1] does not exist 
in the new perturbation theory. 

The calculations of Lifshitz [2] and Lifshitz and Khalatnikov [3] (see also Landau and Lifshitz [4], Chapter 14) are 
redone, and new results of the literature, which were not known at the time when Lifshitz and Khalatnikov developed 
their theory, are used. Combining these results with new insights described in this article, a perturbation theory for 
the open flat or closed Friedmann-Lemaitre-Robertson- Walker (flrw) universe is developed with the properties that 
both the evolution equations and their solutions are invariant under general infinitesimal coordinate transformations 
x'^ x^—S^^{x'^). The Lifshitz-Khalatnikov cosmological perturbation theory has thus been rewritten into a manifestly 
covariant and gauge-invariant perturbation theory, see equations (43). The upshot of the new approach is a possible 
explanation of primeval stars, the so-called (hypothetical) population III stars, independent of the existence of Cold 
Dark Matter (cdm). 

A. Former Insights 

Bardeen was the first to point out the importance of gauge-invariant variables in the construction of a perturbation 
theory. In his seminal article [5] Bardeen demonstrated that the use of gauge-invariant variables ensures that a 
perturbation theory is free of spurious solutions. The article of Bardeen has inspired the pioneering works of Ellis 
et al. [6-8] and Mukhanov et al. [9, 10]. These researchers proposed alternative perturbation theories. 

Although Lifshitz and Khalatnikov were aware of the fact that the system of perturbation equations can be di- 
vided into three independent systems of equations, namely for tensor perturbations (gravitational waves), vector 
perturbations (vortexes), and scalar perturbations, they did not explicitly use the decompositions (17) and (20) found 
by York [11], Stewart [12] and Stewart and Walker [13]. These decompositions make the computations much more 
tractable. 

An important result [11-13] is that the perturbed metric tensor for scalar perturbations can be written in terms of 
two potentials (19). This facilitates the derivation of the Newtonian results (66) and (67) in the non-relativistic limit. 

B. New Insights 

In addition to the results found in the literature, some new insights, described in this article, are used. 

Firstly, barotropic equations of state for the pressure p = p{e), where e is the energy density, are commonly used in 
cosmology [14]. In particular, the linear barotropic equation of state p = we, with w = ^ in the radiation-dominated 
era and w = in the era after decoupling of matter and radiation, describes very well the global evolution of an 
wraperturbcd flrw universe. As will become clear in Section VI, p = we gives an incomplete description of the 
evolution of density perturbations in the perturbed universe. Therefore, realistic equations of state are needed. From 
thermodynamics it is known that both the energy density e and the pressure p depend on the independent quantities 
n and T, i.e., 

£ = e(n,T), p^p{n,T), (1) 

where n is the particle number density and T the temperature. Since T can, in principle, be eliminated from these 
equations of state, a computationally more convenient equation of state for the pressure will be used, namely 

p^p{n,e). (2) 

This, non-barotropic, equation of state has been included in the new perturbation theory (43). Therefore, this theory 
will be referred to as the generalised Lifshitz-Khalatnikov theory. 

Secondly, for scalar perturbations the three perturbed momentum constraint equations (10b) can be rewritten as 
one evolution equation (28b) for the local perturbation (15), to the global spatial curvature R^o)^ (8). Since, in 
first-order, only the irrotational part of the spatial part of the fluid four-velocity is coupled to scalar perturbations, 
the three momentum conservation laws (lOe) can be recast in one evolution equation (28d) for the divergence 
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(16), of the spatial part of the fluid four-velocity. Using the variables and i9(i), the evolution equations for 
scalar perturbations can be written as an initial value problem (28) with four ordinary differential equations and one 
algebraic equation, to be obeyed by the initial values. The only gradient term in this system occurs in equation (28d). 
This is a distinct advantage over former perturbation theories [1, 5-10] which contain metric gradients. 

Thirdly, since the sets of equations (7) and (28) are written with respect to the same system of reference, one can 
draw the conclusion that only three independent scalars, namely the energy density e, the particle number density n 
and the expansion scalar 9, play a role in a perturbation theory. This reduces the number of possible gauge-invariant 
variables considerably, and, thereupon, leads to the insight that there exist two, and only two, unique gauge-invariant 
quantities ef'j and n^',, defined by (39a). 

Finally, in order to show that efl^ and nfl^ are indeed the perturbations to the energy density and particle number 
density, respectively, one has to investigate the basic equations (28) for scalar perturbations together with the defini- 
tions (39a) in the so-called non-relativistic limit. In this limit, which will be defined in exact terms in Section IIIH, 
the set (28) combined with the definitions (39a) reduce to the Newtonian Theory of Gravity, i.e., the Newtonian 
results (66) and (67) show up. Therefore, ef', and nf^) are the real, physical, energy density perturbation and the 
particle number density perturbation, respectively. Consequently, the generalised Lifshitz-Khalatnikov theory (43), 
which is derived from (28) using (39a), is a manifestly covariant and gauge- invariant perturbation theory and describes 
unequivocally the evolution of density perturbations in flrw universes. 

II. RESULTS 

Firstly, it follows from the generalised Lifshitz-Khalatnikov theory that perturbations in the total energy density 
are gravitationally coupled to perturbations in the particle number density (43b), irrespective of the nature of the 
particles (i.e., ordinary matter or CDm) and independent of the scale of the perturbations. 

Secondly, in the radiation-dominated era of a fiat flrw universe small-scale perturbations oscillate with an increas- 
ing amplitude according to (76), whereas the standard evolution equation (104) predicts oscillatory perturbations 
with a constant amplitude. This difference is due to the fact that in the new perturbation theory the divergence 
(16), of the spatial part of the fiuid four-velocity has been taken into account, whereas t9(i) = in the standard theory, 
compare the Newtonian equation (104) with the exact relativistic equations (107). In the evolution of large-scale per- 
turbations 'd(i) plays only a minor role, so that the results of the standard theory are recovered (75). Heat exchange 
does not play a role in the radiation-dominated era. 

Finally, in the era after decoupling of matter and radiation, heat loss. Section III G, of a perturbation, in addition to 
gravity, has a more or less favourable effect on its growth rate, depending on the scale. For large-scale perturbations 
heat loss is unimportant and gravity is the main cause of contraction. Therefore, the generalised Lifshitz-Khalatnikov 
theory corroborates the results for large-scale perturbations (88) of former perturbation theories which did not take 
heat loss into account. Small-scale perturbations benefit more from heat loss during their evolution than large-scale 
perturbations do. It is shown in Section VC that for perturbations with scales around 6.2 pc gravity and heat loss 
combine optimally. Figure 1, resulting in fast growing density perturbations. This may shed new light on the evolution 
of small-scale inhomogeneities in the universe. 

III. GENERALISED LIFSHITZ-KHALATNIKOV PERTURBATION THEORY 

In this section the generalised Lifshitz-Khalatnikov theory of cosmological density perturbations is derived for the 
open, flat or closed flrw universe. To that end, a suitable system of reference must be chosen. 

Due to the general covariance of the Einstein equations and conservation laws, Einstein's gravitational theory is 
invariant under a general coordinate transformation — > x^{x^), implying that preferred coordinate systems do not 
exist (Weinberg [15], Appendix B). In particular, the linearised Einstein equations and conservation laws are invariant 
under a general linear coordinate transformation 

-> a;° - V(t,ic), x' ^ x' - C{t,x), (3) 

where ip(t,x) and ^*(t, a;) are four arbitrary, first-order (infinitesimal) functions of the time {x^ = ct) and space 
[x = {x^,x^,x'^)\ coordinates. Since preferred systems of reference do not exist and since the final result (43) is 
manifestly covariant and gauge- invariant, one may use any suitable and convenient coordinate system to perform the 
calculations. In order to put an accurate interpretation on the new gauge-invariant quantities (39) one needs the non- 
relativistic limit. In the Newtonian Theory of Gravity all possible systems of reference are synchronous. Therefore, 
synchronous coordinates [2-4], i.e., Qqi = 0, will be used in the background as well as in the perturbed universe. For 
convenience, synchronous coordinates are chosen such that coordinate time is equal to proper time, i.e., 500 = 1- 
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A. Basic Equations in Synchronous Coordinates 



In synchronous coordinates the metric of flrw universes has the form 

goo = 1, goi = 0, gtj = -a'^{t)gij{x), (4) 

where a(t) is the scale factor of the universe, and gij is the metric of the three-dimensional maximally symmetric 
subspaces of constant time. The functions ^ and i^' of the general infinitesimal transformation (3) reduce to 

i^^^Pix), f r^+X^N, (5) 

J a-^lT) 

if only transformations between synchronous coordinates are allowed. In (5), ipix) and x*(^) ^^'^ four arbitrary, 
infinitesimal functions of the spatial coordinates. 



1. Background Equations 



The complete set of zeroth-order Einstein equations and conservation laws for an open, flat or closed flrw universe 
filled with a perfect fluid with energy-momentum tensor 

T^^- ^{e+p)u^^u'' ^pg^"-, v = v{n,e), (6) 

is, in synchronous coordinates, given by 

iH^ = ^R(o) + ne^oi + = SttGn/c*, (7a) 

^(0) ~ — 2i/_R(o), (7b) 

£(0) = -3He(o)(l + w), w = p^o)/s^o), (7c) 

^(0) = 0, (7d) 

71(0) = -3Hn(o). (7e) 

The Goi constraint equations and the Gij, i ^ j, dynamical equations are identically satisfied. The Gu dynamical 
equations are equivalent to the time-derivative of the Gqq constraint equation (7a). Therefore, the Gij dynamical 
equations need not be taken into account. In equations (7) A is the cosmological constant, Gn the gravitational 
constant of the Newtonian Theory of Gravity, c the speed of light and w a shorthand notation for the quotient 
P(o)/£(o)- An over-dot denotes differentiation with respect to ct and the sub-index (0) refers to the background, i.e., 
unperturbed, quantities. Furthermore, H = a/a is the Hubble function which is equal to H = |0(o), where 9^o) is the 
background value of the expansion scalar 9 = u'^-f^ with = c~^U'^ the four-velocity, normalised to unity {u^^u^ = 1). 
A semicolon denotes covariant differentiation with respect to the background metric g(o)p,i'- The spatial part of the 
background Ricci curvature tensor Rl„fj and its trace R^o) are given by 

2K 6K 

^(0)7= T^^V R(o) = 5"' K = -l,0,+l, (8) 

where i?(o) is the global spatial curvature. The quantity is the three-divergence of the spatial part of the four- 
velocity ufo). For an isotropically expanding universe the four-velocity is u^g^ = S^o, so that = 0. 

From the system (7) one may infer that the evolution of an unperturbed flrw universe is determined by exactly 
three independent scalars, namely 

e^T^'u^u,, n = mu^, = u''.^, (9) 



where iV^ = nu^ is the cosmological particle current four-vector, which satisfies the particle number conservation law 
iV.^ = 0, (7e), see Weinberg [15], Appendix B. As will become clear in Section HID, the quantities (9) and their 
first-order counterparts play a key role in the evolution of cosmological density perturbations. 
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2. Perturbation Equations 

The complete set of first-order Einstein equations and conservation laws for the open, flat or closed flrw universe 
is, in synchronous coordinates, given by 

ij/i'^'fc + = (fOa) 

h^k\i - h^i\k = 2K(e(o) +P(o))U(i)i, (10b) 

K'j + iHh'j + 5'jHh''k + '2.R\,yj - (10c) 

£(1) + 3i/(e(i) + (£(0) +P(o))6'(i) = 0, (lOd) 

+^<«))"(i) ~ 9\o)Pw\k + 5-ff(£(o) +Pm)u\^^ = 0, (fOe) 
+ 3i?n(i) + n(o)^(i) = 0, (lOf) 

where /i^^ = —g^i^f^i, with /ipo = 0, Hqi = is the perturbed metric, h'^j — g^^^^hkj, and g\l-^ is the unperturbed 
background metric (4) for an open, flat or closed flrw universe. Quantities with a sub-index (1) are the first-order 
counterparts of the background quantities with a sub-index (0) . A vertical bar denotes covariant differentiation with 
respect to 5(o)ij- 

The first-order perturbation to the pressure is given by the perturbed equation of state 

=Pn"(i) +Pe£(i), = (1^) ' ' 

where p„(n,£) and Pe{n,e) are the partial derivatives of the equation of state p(n,e). 
The first-order perturbation to the spatial part of the Ricci tensor (8) reads 

= (3'''^fej)(i) = 9\o)R{i)kj + ^R{o)h^j- (12) 
Using Lifshitz' formula [see Weinberg [f6], equation (10.9.1)] 



rfi)y - ^hSmihiiy + - hij\i), (13) 
and the contracted Palatini identities [see Weinberg [16], equation (10.9.2)] 

-^(i)'y — ^^i)ij\k ~ rfi)ife|j' (14) 

one finds, using g\l)h!^i\j\k — g\i)h^i\k\j, for the trace of (12) 

^(1) = 9\i)ih^k\i\] - h^i\k\j) + \R[o)h*'k- (15) 



Expression (15) is the local perturbation to the global spatial curvature i?(o) due to a local density perturbation. 

Using that u'^g 



Finally, ^(i) is the first-order perturbation to the expansion scalar 9 = u^-,^- Using that u'^„^ = 5^q, one gets 



= ~ i5a)^<,|fe, (16) 

where is the divergence of the spatial part of the perturbed four-velocity u'^iy The quantities (15) and (16) play 
an important role in the derivation of the manifestly covariant and gauge- invariant perturbation theory (43). 



B. Decomposition of the Metric and the Spatial Part of the Fluid Four- Velocity 

Since the work of Lifshitz and Khalatnikov, it is known that for flrw universes the set (10) can be broken down 
into three independent systems of equations, namely a system for gravitational waves (tensor perturbations), for 
vortexes (vector perturbations), and for scalar perturbations. In addition they showed that only scalar perturbations 
are coupled to density perturbations. Lifshitz and Khalatnikov used spherical harmonics to classify the three types 
of perturbations. A disadvantage of their approach is that using spherical harmonics is computationally demanding. 
A new insight, gained in 1974 by York [11], and in 1990 elucidated by Stewart [12], makes the computations much 
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more tractable. York and Stewart showed that the symmetric perturbation tensor hij can uniquely be decomposed 
into three parts, i.e., 

h^^^hl^+h^^+hl^, (17) 

where the scalar, vector and tensor perturbations are denoted by ||, _L and *, respectively. The constituents have the 
properties 

hi^^o, h':^ = o, /it-|, = o. (18) 

Moreover, York and Stewart demonstrated that the components /ly^- can be written in terms of two independent 
potentials (t>{t,x) and C,{t,x), namely 

% = |(0<5^+Cl^,). (19) 

Finally, Stewart also proved that the spatial part of the perturbed four-velocity can uniquely be decomposed into 
two parts 

■"(1) = "(1)11 +M(i)-L, (20) 

where the constituents have the properties 

V • tt,i) = V • tt(i)|| , V X tt,i) = V X M(i)_L, (21) 

with V* = g^-'dj, the generalised vector differential operator. 

The three different kinds of perturbations will now be considered according to the decompositions (17) and (20), 
and their properties (18) and (21). 

For tensor perturbations the properties (18) imply that -R(i)* — 0, as follows from (f5). Using this and (18), 
equations (lOa)-(lOc) imply that tensor perturbations are not coupled to £(1), and it(i). 

The perturbed Ricci tensor, R(i-,ij, being a symmetric tensor, should obey the decomposition (17) with the properties 
(18), namely -Rj^jj^^, = 0. This implies that, see (15), h^^ must obey h^j^^^ — 0, in addition to /i^^j. = 0. From (lOa) 
and the trace of (10c) it follows that vector perturbations are not coupled to £(1) and p^iy Raising the index i of 
equation (10b) with g'^l^ one finds that this equation reads for vector perturbations 

/i^,^ + 2i//ig^ = 2k(£(o, +P(o,X,, (22) 

where it is used that i?^^) = —2Hg\ly Taking the covariant derivative of (22) with respect to the index j one finds with 
the additional property h*^^m = that equation (22) reduces to u'^^^^j = 0, implying with (21) that the constituent 
M(i)_L is coupled to vector perturbations. 

Since both 7^ and -R(i)|| 7^ 0, scalar perturbations are coupled to £(i) andp(i). It will now be demonstrated that 
is coupled to scalar perturbations, by showing that equation (10b) requires that the rotation of vanishes, 
if the metric is of the form (19). Differentiating (10b) covariantly with respect to the index j and substituting (19) 
yields 

2'/>|i|j +C}^\k\t\j -i}^\i\k\3 = Kc^(e(o) +P^o))u^^)^\J. (23) 
Interchanging i and j and subtracting the result from (23) one gets 

C}''\i\k\] - C^^\j\k\i = -Kc^(£(o) +P(o))("(i)i|i - W(i)jk)- (24) 
By rearranging the covariant derivatives, (24) can be cast in the form 

(C'^Kifcij - C'^Kiiifc) - {Q'^UMi ~ C'^'bUifc) 

+ (C'''|i|j - C'''bK)|fe = -Kc2(£(0) +P(o))(u(i)i|j - (25) 
Using the expressions for the commutator of second order covariant derivatives (Weinberg [16], Chapter 6, Section 5) 

4* A'i J^k Ak ryi f>i p« f?^ ^9^^"^ 

^ j\p\q ^ j\g\p ~ ^ fc^(Q)jpg ^ j^mkpqi ^ \p\q ^ \q\p ~ ^ ^{a)kpq^ x'^^) 
and substituting the background Riemann tensor for the three-spaces of constant time 

- K {5\~g,i - 5\~g,k) , K ^ -I, 0, +1, (27) 

one finds that the left-hand side of equation (25) vanishes identically, implying that the rotation of u^^ is zero. 
Therefore, only is coupled to scalar perturbations. 
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C. First-order Equations for Scalar Perturbations 

Since scalar perturbations, i.e., perturbations in £(i) and are only coupled to k'^ and one may replace 

in (10)-(16) h'j by hy and by to obtain perturbation equations which exclusively describe the evolution 

of scalar perturbations. From now on, only scalar perturbations are considered, and the subscript |j will be omitted. 
Using the decompositions (17) and (20) and the properties (18) and (21), one can rewrite the evolution equations for 
scalar perturbations in the form 

2H{9^,, - 79(1,) - = Ke(i„ (28a) 
R^,, + - 2Ke(o,(l + w),9,,, + lR^o,iO,^, - ^^^,) = 0, (28b) 

£(1) + 3i?(£(i) +P(i,) + £(0,(1 + w)6'(i) = 0, (28c) 

+ Hi2 - Sm.^ + - (28ci) 

n(i) + 3iJn(i) + n(o)^^(l) = 0. (28e) 

The set (28), which is the perturbed counterpart of the set (7), consists of one algebraic equation (28a) and four 
ordinary differential equations (28b)-(28e) for the five unknown quantities £(i), 71(1, , i?^) and 9^iy The quantity 
/3(i) is defined by = p^o■,/e^oy Using that p^a) — Pnn^o) + Pe^io) and the conservation laws (7c) and (7e) one gets 

£(0)(l + w) 

Finally, the symbol denotes the generalised Laplace operator with respect to the three-space metric gij, defined 
byV2/ = 5'V|.b-- 

The derivation of the basic equations (28) for scalar perturbations will now be given. Eliminating h'^k from (10a) 
with the help of (16) yields the algebraic equation (28a). 

Multiplying both sides of equation (10b) by glgj and taking the covariant derivative with respect to the index j, 
one finds 

9\i){h^k\i\j - h^i\k\j) = 2k(£(o) +P(o))i9(i), (30) 

where also (16) has been used. The left-hand side of (30) will turn up as a part of the time-derivative of the curvature 
_R(i). In fact, differentiating (15) with respect to time and recalling that the background connection coefhcients rj^^^j^ 

are for flrw metrics (4) independent of time, one gets, using also g\l-^ = —2Hg\l-^ and (7b), 

= -"^HRfj^^ + g\i){^^ k\i\j — h^i\k\j) + \R{0)h^k- (31) 

Combining (30) and (31) and using (16) to eliminate h^k yields (28b). Thus, the G°j^jj momentum constraint equations 
(10b) have been recast into one equation (28b) for the local spatial curvature due to a density perturbation. 

For i 7^ j equations (10c) need not be considered, since they are not coupled to scalar perturbations. Taking the 
trace of (10c) and eliminating the quantity h!^k with the help of (16), one arrives at 

^(1, - 79,1, + 6ff ((?(!, - 79,1,) - = §^£(£(1) - P(i)). (32) 
Using (28a) to eliminate the second term of (32) yields for the trace of (10c) 

0,1, - 79,1, + \Rii) = -i'«(e(i) +P(i))- (33) 
This equation is equal to the time-derivative of the constraint equation (28a), which reads 

2ij(0,i, - 79,1,) + 2i7(^,i, - 79(1,) - i7?(i, - K£(i,. (34) 

Eliminating the time-derivatives -ff, 7?(i, and £,i, with the help of (7a)-(7c), (28b) and (28c), respectively, yields 
the dynamical equation (33). Consequently, for scalar perturbations the dynamical equations (10c) need not be 
considered. 

Finally, taking the covariant derivative of (10c) with respect to the metric (?,o)ij and using (16), one gets 

(£(0) +P(o))i9(i) -ff(o,P(i)|/cK + 5i?(£,o)+P(o))^?(i) =0, (35) 
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where it is used that the operations of taking the time-derivative and the covariant derivative commute, since the 
background connection coefficients T^^^^^j are independent of time for flrw metrics. With (4), V^/ = g^-' f\i\j and (7c) 
one can rewrite (35) in the form 

+ i? f 2 - 3^) + ^ = 0, (36) 

\ ^(0) / £(0) + P(o) a 

Using the definitions w = P(o)/£(o) and = p^o)/£(o) one arrives at equation (28d). 

This concludes the derivation of the system (28). As follows from its derivation, this system is, for scalar perturba- 
tions, equivalent to the full set of first-order Einstein equations and conservation laws (10). 



D. Unique Gauge-invariant Cosmological Density Perturbations 

The background equations (7) and the perturbation equations (28) are both written with respect to the same system 
of reference. Therefore, these two sets can be combined to describe the evolution of the five background quantities 
6'(o) = 37J, i?(o), £(0), ^^(0) = 0, 71(0), and their first-order counterparts -R(i), £(1), ■(?(!), 7t,(i). Just as in the background 
case, one again comes across the three independent scalars (9). Consequently, the evolution of cosmological density 
perturbations is described by the three independent scalars (9) . A complicating factor is that the first-order quantities 
£(1) and 71(1), which are supposed to describe the energy density and the particle number density perturbations, have 
no physical significance, as will now be established. 

A first-order perturbation to one of the scalars (9) transforms under a general (not necessarily between synchronous 
coordinates) infinitesimal coordinate transformation (3) as 

S^,)it,x) ^ S^,y{t,x) +ij{t,x)S,o^{t), (37) 

where S^o) and Sf^ are the background and first-order perturbation of one of the three scalars S = e,n,9. In (37) 
the term S = 'tpS^o) is the so-called gauge mode. The complete set of gauge modes for the system of equations (28) is 
given by 



£(i)=V'£(o), "(1) = '^"■(o), = V'6'(o), (38a) 

(38b) 



„2 2^mV 



where ip = ipix) in synchronous coordinates, see (5). The quantities (38) are mere coordinate artifacts, which have no 
physical meaning, since the gauge function ip^x) is an arbitrary (infinitesimal) function. Equations (28) are invariant 
under coordinate transformations (3) with (5), i.e., the gauge modes (38) are solutions of the set (28). This property 
combined with the linearity of the perturbation equations, implies that a solution set (e(i), n(i), 0(i), i9(i), i?(i)) can 
be augmented with the corresponding gauge modes (38) to obtain a new solution set. Therefore, the solution set 
(£(1), 71(1), 0(1), ■(^(i), i?(i)) has no physical significance, since the general solution of the set (28) can be modified by an 
infinitesimal coordinate transformation. This is the notorious gauge problem of cosmology. 

In this article, the cosmological gauge problem has been resolved. To that end, the perturbation equations (10) have 
first been rewritten into the form (28) in order to isolate the scalar perturbations from the vortexes and gravitational 
waves. This reduces the number of possible gauge-invariant variables substantially, since one needs to consider only 
the three independent scalars (9). All scalar perturbations transform under the general infinitesimal transformation 
(3) in the same way according to (37). This implies that one can combine two independent scalars to eliminate 
the gauge function ij}{t,x). With the three independent scalars (9), one can make (2) = 3 different sets of three 
gauge-invariant variables. In each of these sets exactly one gauge-invariant quantity vanishes. As will be shown in 
Section IIIH, the only set for which the corresponding perturbation theory yields the Newtonian results (66) and (67) 
in the non-relativistic limit is given by 



77o 



£(1) ~ £(1) /, (1)' (1) ~ ' (1) n 



qgi _ n _ "(0) 
7(1) — ^(1) ■ ■^(1) 



1) = 0. (39b) 
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It follows from the general transformation rule (37) that the quantities (39) are invariant under the general infinitesimal 
transformation (3), i.e., they are gauge-invariant, hence the superscript 'gi'. The definitions (39) imply that the gauge- 
invariant counterpart Bfl-^ of the gauge dependent variable 6'(i) ^ vanishes automatically. The physical interpretation 
of (39b) is that, in first-order, the global expansion 0(o) — 3H is not affected by a local density perturbation. 

The quantities (39a) have two essential properties, which gauge-invariant quantities used in former perturbation 
theories [5-10] do not have. Firstly, due to the quotient S^o)/(^(o) of time derivatives, the definitions (39) are independent 
of the definition of time. As a consequence, the evolution of ef', and n^', will only be determined by their propagation 
equations. Secondly, the quantities (39a) do not contain spatial derivatives, so that unnecessary gradients [1] do not 
occur in the final equations (43). 

The quantities (39a) are completely determined by the background equations (7) and their perturbed counterparts 
(28). In principle, these two sets can be used to study the evolution of density perturbations in FLRW universes. 
However, the set (28) is still too complicated, since it also admits the non-physical solutions (38). The aim will 
be a system of equations for e^', and nfi, that do not have the gauge modes (38) as solution. In other words, a 
perturbation theory will be constructed for which not only the differential equations are invariant under general 
infinitesimal coordinate transformations (3), but also their solutions. Such a theory will be referred to as a manifestly 
covariant and gauge-invariant perturbation theory. The construction of this theory will be the subject of the next 
subsection. 



E. Manifestly Covariant and Gauge-invariant Perturbation Theory 

In this section the derivation of the new perturbation theory is given. Firstly, it is observed that the gauge dependent 
variable 9^l^ is not needed in the calculations, since its gauge-invariant counterpart 9fl^, (39b), vanishes identically. 
Eliminating 6'(i) from the differential equations (28b)"(28c) with the help of the (algebraic) constraint equation (28a) 
yields the set of four first-order ordinary differential equations 



R 

+ 2i/i?(i) - 2Ke(o)(l + w)^9(i) + (k£(i, + ii?,!,) = 0, 

£(1) + 3iJ(e(i) + P(i,) + £(0) (1 + W) + ^ (k£(i) + ii?(i)) 



= 0, 



£(0)(l + w) a2 



0, 



= 0, 



(40a) 
(40b) 

(40c) 

(40d) 



for the four quantities £(i), nji), and 

Using the background equations (7) to eliminate all time-derivatives and the first-order constraint equation (28a) 
to eliminate the gauge-invariant quantities (39a) become 



£(i)i?(0) ~ 3£(o)(l + w){2Hd,,^ + ifl(i)) 
i?(0) -I- 3k£(o)(1 + w) 



'(1) — "(l) 



3n(o)(K£(i) -\- 2iJ?9(i) ^i?(i)) 

i?(0) -I- 3k£(o)(1 + w) 



(41a) 
(41b) 



These quantities are completely determined by the background equations (7) and the first-order equations (40). In 
the study of the evolution of density perturbations, it is convenient not to use efl^ and nfl^ directly, but instead their 
corresponding contrast functions and i5„ 

5,{t,x) ^ S^{t,x) ^ (42) 

The system of equations (40) for the four independent quantities £(i), n(i), 'd^i^ and will now be rewritten, with 
the help of a computer algebra system, into a new system of equations for the two independent quantities and 
In this procedure it is explicitly assumed that p ^ 0, i.e., the pressure does not vanish identically. The case p — >■ 
will be considered in Section IIIH on the non-relativistic limit. The final result is the generalised manifestly covariant 
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and gauge-invariant perturbation theory for the open, flat or closed FLRW universe 



1 d 

cdt 



Sn 



blSe + b2Se 



1 



1 



e(o)(l + w) 



1 



These are two linear differential equations for the two unknown quantities and $„ 
and equation (7c) that the time-derivative of w = p^o)/ ^(o) is 



w 



(43a) 
(43b) 

It follows from = P(o)/^(o) 

(44) 



Defining p„ 
by 



d'^p/d'n? and p^n = d'^p/dedn and using (44) the coefficients 6i, 62 and 63 of equation (43a) are given 



b2 



K£(o)(l + W) _ P_ 



1 

3H 



2i7(l + 3/32 



i?(o) + 3Ke(o)(l + w) 



iK£(o)(l + w){l + 3w) + H^{l-3w + 6/3^(2 + 3w)) + 6H^ (^i 



K£(0) (1 



i?(o) + 3k£(o)(1 + w) 



63 =' 



R 

'"' + 3Ke(„)(l -I- ly) 

R^o) + 3k£(o)(1 + w) 



T R(t 



e(0)Pen(l + W) + +Pn{Pe 



P^) + n^o)Pnn 



Pn 



2 -^(0) 



(45a) 



(45b) 



(45c) 



The equations (43) have been checked (see attached maple file) using a computer algebra system, as follows. Substi- 
tuting the contrast functions (42) in equations (43), where £f') and nfl^ are given by (41), and subsequently eliminating 
the time-derivatives of £(o), n^gy, H, i?(o) and £(i), n(i), with the help of equations (7) and (40), respectively, 

yields two identities for each of the two equations (43). 

The equations (43) are two differential equations for the two independent and gauge-invariant quantities 5^ and 
5n- It follows from equation (43b) that perturbations in the total energy density are gravitationally coupled to 
perturbations in the particle number density if p„ = {dp/dn)^ < 0, or, equivalently, Pe = {dp/de)n > /3^, see (29). 
This is the case in a flrw universe in the radiation-dominated era and in the era after decoupling of matter and 
radiation. This coupling is independent of the nature of the particles, i.e., it holds true for ordinary matter as well as 
CDM. 

The system of equations (43) is equivalent to a system of three first-order differential equations, whereas the original 
set (40) is a fourth-order system. This difference is due to the fact that the gauge modes (38), which are solutions of 
the set (40), are completely removed from the solution set of (43): one degree of freedom, namely the gauge function 
■0, has disappeared altogether. 

The background equations (7) and the new perturbation equations (43) constitute a set of equations which enables 
one to study the evolution of small fluctuations in the energy density (5^ and the particle number density (5„ in an 
open, flat or closed flrw universe with A ^ and filled with a perfect fluid with a non-barotropic equation of state 
for the pressure p = p{n, s). 



F. Gauge-invariant Pressure and Temperature Perturbations 



The gauge-invariant pressure and temperature perturbations, which are needed in the forthcoming sections, will 
now be derived. 

From the equation of state (2) for the pressure p — p(n, e) it follows that 

P(o)=Pnn^o)+Pe£m, ^"=(1^) ' ^'^ (If ) ■ ^^^^ 

Multiplying both sides of this expression by ^(ij/^fo) and subtracting the result from given by (11), one gets, using 
also (39a), 
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Hence, the quantity defined by 



— P(i) A ''(1)1 



(48) 



is the gauge-invariant pressure perturbation. 

From the equation of state (1) for the energy density e = e{n,T) it fohows that 



de 



de_\ 



de 



de_ 
df 



(49) 



Multiplying £(0, by 6l(i,/0 (0) and subtracting the result from one finds, using (39a), 



de 
df 



rp ^(0) a 
A ''(1) 



(50) 



implying that the quantity defined by 



T 

rrigl — rp _ ^(0) a 

— A I'd)! 



(51) 



is the gauge-invariant temperature perturbation. The expressions (48) and (51) are both of the form (39). 



G. Diabatic Density Perturbations 

In this section equations (43) will be linked to thermodynamics and it will be shown that, in general, density 
perturbations evolve diabatically, i.e., they exchange heat with their environment during their evolution. 
The combined First and Second Law of Thermodynamics is given by 

dE ^ TdS - pdV + ^idN, (52) 

where E, S and N are the energy, the entropy and the number of particles of a system with volume V and pressure p, 
and where fi, the thermal — or chemical — potential, is the energy needed to add one particle to the system. In terms 
of the particle number density n = N/V, the energy per particle E/N = e/n and the entropy per particle s = S/N 
the law (52) can be rewritten 

d (^iv) = Td{sN) - pd (^^^ + fidN, (53) 

where e is the energy density. The system is extensive, i.e., S{aE, aV, aN) — aS{E, V, N), implying that the entropy 
of the gas is S = {E + pV — ^N)/T. Dividing this relation by N one gets the so-called Euler relation 

M = - Ts. (54) 

n 

Eliminating /x in (53) with the help of (54), one finds that the combined First and Second Law of Thermodynamics 
(52) can be cast in a form without /x and N , i.e., 

Tds = df-) +pdf-). (55) 



From the background equations (7) it follows that S(o) = 0, implying with (37) that S(i) — sfl-^ is automatically 
gauge- invariant . 

The thermodynamic relation (55) can, using (42), be rewritten in the form 



gi _ £(o)(l + w) 

J(0)A(i) — 



5n 



l + 'W_ 

Thus, the right-hand side of (43a) is related to local perturbations in the entropy. 



(56) 
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Adiabatic perturbations do not exchange heat with their environment, so that T^o)sfl^ = 0. This impUes with (56) 
that {l + w)dn — 5e = 0. Muhiplying this expression by 3He^o)n^o) and substituting (42) one finds from the background 



conservation laws (7c) and (7e) that the adiabatic condition sfl^ 
the latter expression becomes 



de_ 
df 



't(0)J(l) 



reads 



(-(l) J- (0) 



-£(0)71^1) — 0. Using that e — e{n, T) 



(57) 



Since n and T are independent quantities and since in a non-static universe one has njo) 7^ and r(o) 7^ 0, the adiabatic 
condition (57) is satisfied if, and only if. 



df 



= 0, 



(58) 



implying that e = e{n). In particular, in the non-relativistic limit, where e = nmc^ and p = 0, density perturbations 
are adiabatic. Hence, in all other cases where p — p{n^ e) local density perturbations evolve diabatically. 



H. Non-relativistic Limit 

In Section III D it has been shown that the two gauge-invariant quantities efl^ and nfl^ are unique. It will be 
demonstrated that in the non-relativistic limit equations (28) combined with (39) reduce to the results (66) and 
(67) of the Newtonian Theory of Gravity and that the quantities efl^ and nfl^ become equal to their Newtonian 
counterparts. 

The non-relativistic limit is defined by three requirements, Carroll [17]. Firstly, the gravitational field should be 
weak, i.e., can be considered as a perturbation of flat space. Secondly, the particles are moving slowly with respect 
to the speed of light. Finally, the gravitational field of a density perturbation should be static, i.e., it does not 
change with time. This definition of the non-relativistic limit, which is essential to put an accurate interpretation on 
the quantities efl^ and nfl^, has not been used in former perturbation theories [5-10] to explain the meaning of the 
gauge-invariant variables. 

In a first-order cosmological perturbation theory the gravitational field is already weak. In order to meet the first 
requirement, a flat (i?(o) — 0) flrw universe is considered. Using (19), the local perturbation to the spatial curvature 
(15) reduces for a flat flrw universe to 



4 



(59) 



where is the usual Laplace operator. Substituting this expression into the perturbation equations (28), one gets 



47rGN 



1 



47rGN 

^4 



£(0)(1 -I- U')t?(i) = 0, 



£(1) + 3i7(£(i, +P(i)) + £(o)(l + w)0(i) = 0, 

1 V2p(„ 



+ i/(2 - 3/3^)z9(i) 
ri(i) -I- 3Hn,r) + n^a)0^l 



£(o)(l 
= 0, 



0, 



(60a) 

(60b) 

(60c) 

(60d) 
(60e) 



where (39a) has been used to eliminate £(i) from the constraint equation (28a). 

Next, the second requirement will be implemented. Since the spatial part of the fluid four- velocity is gauge 
dependent with a physical component and a non-physical gauge part, the second requirement should be defined by [18] 



(1) physical 



r^^TP — v n 

^ "-^ (1) physical ' "i 



(61) 



i.e., the physical part of the spatial part of the fiuid four- velocity is negligible with respect to the speed of light. In 
this limit, the mean kinetic energy per particle |m(t)^) = "^k-oT — > is very small compared to the rest energy m(? 
per particle. This implies that the pressure p = nk^T — )■ (n 7^ 0) is vanishingly small with respect to the rest 
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energy density nmc^ . Substituting p = into the momentum conservation law (lOe) yields, using also the background 
equation (7c) with w = P(o)/^(o) — ^ 0, 

v{,^=-2Hv^,^. (62) 

Since the physical part of w^j^j vanishes in the non-relativistic limit, the general solution of equation (62) is exactly 
equal to the gauge mode (38b) 

u\,,{t,x) = --^^g'Hx)dk^{x), (63) 

where it is used that H = a/a. Thus, in the limit (61) one is left with the gauge mode (63) only. Consequently, one 
may, without losing any physical information, put the gauge mode ttj^j equal to zero, implying that dkip = 0, so that ij^ 
is a constant in the non-relativistic limit. Substituting dk4' = into (5) one finds that the relativistic transformation 
(3) between synchronous coordinates reduces in the limit (61) to the (infinitesimal) transformation 

x' ^ x' -x'{x), (64) 

where is an arbitrary constant and x*(3;) are three arbitrary functions of the spatial coordinates. In the non- 
relativistic limit time and space transformation are decoupled: time coordinates may be shifted and spatial coordinates 
may be chosen arbitrarily. The residual gauge freedom -ip and x'(^) ^^ot come as a surprise, since the Newtonian 
Theory of Gravity is invariant under the gauge transformation (64). 
Substituting — and p — Q into the system (60), one gets 

V^^^^a^ef', (65a) 

V^^ = 0, (65b) 
£(1) -H3i7£(i, -l-e(o)0(i) = 0, (65c) 
n(i, -I- 37Jn(i) -I- »^(o)6'(i) = 0. (65d) 

The constraint equation (65a) can be found by subtracting ^9^i^/H times the time-derivative of the background 
constraint equation (7a) with i?(o) = from the constraint equation (60a) and using that 0(o) = 3i?. Note that the 
cosmological constant A need not be zero. 

Since i?;!, = 0, i.e., no fluid flow, density perturbations do not evolve, so that their gravitational field is static 
(65b). Consequently, a'^{t)efl^{t,x) in (65a) should be replaced by a^(to)£fi) (^oj a^)- Defining the potential (p{x) = 
4>{x) / [to) , equations (65a) and (65b) imply 

VV(a;) = 47rGNP<„(a;), p,„(x) ^ ^M^^ (66) 

which is the Poisson equation of the Newtonian Theory of Gravity. With (66) the third requirement for the non- 
relativistic limit, i.e., a static gravitational field, has been satisfied. The universe is in the non-relativistic limit not 
static, since H ^ and iJ ^ 0, as follows from the background equations (7) with w = and R(o) — 0. In the 
non-relativistic limit a local density perturbation does not follow the global expansion of the universe. 

Equations (65c) and (65d) have no physical significance since the only solutions of these equations are the gauge 
modes (38a) with constant tp, (64). Therefore, the latter two equations are not part of the Newtonian Theory of 
Gravity and need not be considered. The potential C which occurs by (19) in -R(i), (15), and (16), in the general 
relativistic case, drops from the perturbation theory in the non-relativistic limit. Consequently, one is left with one 
potential (p{x) only. 

The expression (41a) reduces in the non-relativistic limit to efl^ = —R^i)/{2k), which is, with (59) and (65b), 
equivalent to the Poisson equation (66). Expression (41b) reduces in the non-relativistic limit to the familiar result 



(67) 



where it is used that, in the non-relativistic limit, e^l^ = n^i^rrK? and £(o) — n^o^m(? . 

It has been shown that equations (28) combined with the unique gauge-invariant quantities (39) reduces in the non- 
relativistic limit to the Newtonian results (66) and (67). Consequently, efi, and nf^^ are the real, physical perturbations 
to the energy density and particle number density, respectively. 
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IV. EXAMPLE: THE FLAT FLRW UNIVERSE 



In this section analytic solutions of equations (43) are derived for a flat (i?(o) = 0) flrw universe with a vanishing 
cosmological constant (A = 0) in its radiation-dominated phase and in the era after decoupling of matter and radiation. 



A. Radiation-dominated Era 



In the radiation-dominated epoch one has e = obT^, where ab is the black body constant and T-y the radiation 
temperature. The pressure is p — |e, so that p„ = 0, pg = |, implying, with (29), that = The perturbation 
equations (43) reduce to 



3 a2 



<5e - H5, 



(68a) 
(68b) 



Equation (68b) implies that the difference (5„ — f^e depends only on the spatial coordinates. Since (68a) has no source 
term, heat exchange is unimportant in the radiation-dominated era. Therefore, one may use the adiabatic condition 
Sn = jSe, see Section IIIG. 

Equation (68a) will now be rewritten in dimensionless quantities. The solutions of the background equations (7) 
are given by 



H(xt-\ 



£(0) (X t 



71(0) oc t a (xt 



1/2 



implying that Tfo^^^ cx a ^. The dimensionless time r is defined by r = t/tQ. Since H = d/a, one finds that 



1 



mitotic 



k = 1,2. 



c^At^ [c<o 

Substituting 5e{t,x) — (5£(i, q) exp(i<7 • x) into equation (68a) and using (70) yields 

" '2 1 - 



4r ^ 2t2 



(5. =0, 



T> 1, 



where a prime denotes differentiation with respect to r. The parameter /Xr is given by 

2tt 1 1 



Air 



Ao if (to) 73' 



Ao = Aa(to), 



(69) 



(70) 



(71) 



(72) 



with Aa(to) the physical scale of a perturbation at time to, and |q| = 27r/A. To solve equation (71), replace r by 
X = i-LrV^- After transforming back to r, one finds 



Se{T, q) = Ai{q) sin {^JL,^/T) + ^2(9) cos {^JL^^/T) 
where the 'constants' of integration Ai{q) and A2{q) are given by 



M{q) = SeitQ,q) sin/ii- - 
^2(9) = Ss{to,q) cos Mr + 



cos /if 
Mr 

sin /ii- 

Mr 



Se{tQ,q) - 



Se{to,q) 
H{ta) 

Seito,q) 
H{to) 



(73) 

(74a) 
(74b) 



For large-scale perturbations (A — > 00), it follows from (73) and (74) that 



Seito) - 



Hito) 



t 



H{to) 



1 

t \ 2 



(75) 
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The energy density contrast has two contributions to the growth rate, one proportional to t and one proportional to 
t^/^. These two solutions have been found, with the exception of the precise factors of proportionality, by a large 
number of authors, see Lifshitz and Khalatnikov [3], (8.11), Adams and Canuto [19], (4.5b), Olson [20], page 329, 
Peebles [21], (86.20), Kolb and Turner [22], (9.121) and Press and Vishniac [23], (33). Consequently, the generalised 
Lifshitz-Khalatnikov theory corroborates for large-scale perturbations the results of the literature. 

A new result is that small-scale perturbations (A — ?■ 0) oscillate with an increasing amplitude according to 



1 

S,{t,q) w (5e(to,q) ( cos 



^h■ - Mr 



(76) 



as follows from (73) and (74). These small-scale perturbations manifest themselves as small temperature fluctuations 
in the cosmic background radiation. 

By virtue of equation (68b) particle number density fluctuations (5„ are coupled to fluctuations in the energy 
density. Since equation (68b) is independent of the nature of the particles and CDM interacts only via gravity with 
ordinary matter and radiation, the fluctuations in CDM are gravitationally coupled to fluctuations in the energy 
density. In other words, 5n — z^e holds true for ordinary matter as well as CDM. Consequently, in the radiation- 
dominated universe CDM does not contract faster than ordinary matter, so that star formation can only commence 
after decoupling. 

The standard equation (104) with w = ^ predicts oscillating density perturbations with a constant amplitude. As 
will be explained in Section VI, this outcome is to be expected since in the standard equation the divergence of the 
spatial part of the fluid four- velocity, 'd(i), has not been taken into account. If iS^^y = there is no fluid flow, implying 
that a density perturbation does not grow. In the new theory (68) 7^ 0, so that oscillating density perturbations 
with an increasing amplitude (76) are found. 



B. Era after Decoupling of Matter and Radiation 

Once protons and electrons combine to yield hydrogen, the radiation pressure becomes negligible, and the equations 
of state (1) become those of a non-relativistic monatomic perfect gas 



e{n, T) = mn(? + ^nk^T, pin, T) = nk-QT^ 



(77) 



where fee is Boltzmann's constant, m the mean particle mass, and T the temperature of the matter. It is assumed 
that the CDM particle mass is larger than or equal to the proton mass, mcDM > "t-h, implying that for the mean 
particle mass m one has mc^ 3> k^T , so that w = P(o}/^(o) ^ 1- Therefore, as follows from the background equations 
(7a) and (7c), one may neglect the pressure nfeeT and the kinetic energy density ^nkBT with respect to the rest-mass 
energy density nmc^ in the unperturbed universe. However, neglecting the pressure in the perturbed universe yields 
density perturbations with a static gravitational field as has been demonstrated in Section IIIH. 

Eliminating T from (77) yields p{n,e) = |(e — nmc^), so that p^ = {dp/de)n = | and p„ = {dp/dn)^ — —^mc^. 
Substituting p„, p^ and (77) into (29) on finds, using mc^ ^ ^bT, 



Vs{t) 



5 fcBT(o,(t) 



^(0) oc a 



(78) 



with Vs the adiabatic speed of sound and T(o) the matter temperature. The fact that T(o) cx a ^ follows from the 
equations of state (77) and the conservation laws (7c) and (7e). This, in turn, implies with (78) that $/P = —H. The 
system (43) can now be rewritten as 



4 + 3HS, 
1 d 



'dt 



2 V 

3 a 



2 {Sn - Se), 



(79a) 
(79b) 



where w <^ 1 and 1 have been neglected with respect to constants of order unity. 

From equation (79b) it follows that 



(5e cx a 



(80) 
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where it is used that H = a/a. 

Using that k-QT^o^ <C mc^, one finds for the perturbed counterparts of (77) 



d„ - de « j-Ot, 

2 mc^ 



<5p = 6n + 5j 



(81) 



where (5p is the relative pressure perturbation defined by 5p = pfl jp^o) and 5t is the relative matter temperature 
perturbation defined by 5t = T^^/T^o), see (48) and (51). Combining (78) and (80) one finds from (81) that 6t is 
nearly constant, i.e., 



STit,x) sa 5Tito,x), 



(82) 



to a very good approximation. 

From (78) and (81) it follows that the source term of equation (79a) is of the same order of magnitude as the term 
with coefficient f3'^. It follows from (45c) that the source term vanishes for barotropic equations of state p = p{e). As 
will become clear in Section V, the source term of equation (79a) is crucial for the understanding of star formation 
in the early universe after decoupling. That is why the realistic equation of state p — p{n, e) has been incorporated 
from the outset in the perturbation theory (43). 

Equation (79a) will now be rewritten in dimcnsionless quantities. The solutions of the background equations (7) 
are given by 



oc t 



cx t 



a (X t 



2/3 



(83) 



where the kinetic energy density and pressure have been neglected with respect to the rest-mass energy density. The 
dimcnsionless time r is defined by r = t/tQ. Using that H = a/a, one gets 



d" 



1 

Ctn 



dr*^ 



[i^(^o)]' 



1,2. 



(84) 



Substituting (?£(*, a;) = ^^(t, q) exp(i<7 • x) and 5„(t,a;) = (5„(t, q) exp(iq • x) into equations (79) and using (81) and 
(84) one finds that equations (79) can be combined into one equation 



5': + -5' + 



4 



9^8/3 



10 

9^2 



4 Mm 
15 t8/3 



r > 1, 



where a prime denotes differentiation with respect to r. The parameter /im is given by 

Ao = Aa(to), 



27r 1 Vs{to) 



(85) 



(86) 



with Aa(io) the physical scale of a perturbation at time to, and |q| = 2tt/X. To solve equation (85) replace t by 
X = 2ii^T~^^^. After transforming back to r, one finds for the general solution of (85) 



S,{T,q)= Bi(qr)J+7(2/iinr" 



-l/3\ 



B2{q)J_7{2fi^T 



-1/3] 



-1/2 



5t2/3 



(87) 



where J±7/2{x) are Bessel functions of the first kind and Bi{q) and B2{q) are the 'constants' of integration. 
In the large-scale limit X ^ oo terms with vanish, so that the general solution of equation (85) is 



Seit) 



H{to) 



Sejto) 

Hito) 



Thus, for large-scale perturbations the initial value ST{to,q) does not play a role during the evolution: large-scale 
perturbations evolve only under the influence of gravity. These perturbations are so large that heat exchange does 
not play a role during their evolution in the linear phase. For perturbations much larger than the Jeans scale (i.e., 
the peak value in Figure 1), gravity alone is insufficient to explain star formation within 13.75 Gyr, since they grow as 
(5e cx i^^'^. The solution proportional to t'^^^ is a standard result. Since ^£ is gauge-invariant, the standard non-physical 
gauge mode proportional to is absent from the new theory. Instead, a physical mode proportional to i"^/^ is 
found. This mode has also been found by Bardeen [5], Table I, and by Mukhanov et al. [9], expression (5.33). In order 
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to arrive at the t~^/^ mode, Bardeen has to use the 'uniform expansion gauge.' In the generahsed Lifshitz-Khalatnikov 
theory the Hubble function is automatically uniform, (39b), without any additional gauge condition. Consequently, 
is in agreement with results given in the literature. 
In the small-scale limit A — 0, one finds 



5,{t,q) 



4STitQ,q) 



UTitQ,q) + Se{to,q) 



2/i„ 



2Mm ( ^ 



(89) 



Thus, density perturbations with scales much smaller than the Jeans scale oscillate with a decaying amplitude which 
is smaller than unity: these perturbations are so small that gravity is insufficient to let perturbations grow. Heat loss 
alone is not enough for the growth of density perturbations. Consequently, perturbations with scales much smaller 
than the Jeans scale will never reach the non-linear regime. 

In the next section it is shown that for density perturbations with scales of the order of the Jeans scale, the action 
of both gravity and heat loss together may result in massive stars several hundred million years after decoupling of 
matter and radiation. 



V. EVOLUTION OF SMALL-SCALE INHOMOGENEITIES AFTER DECOUPLING 



In this section it is demonstrated that the generalised Lifshitz-Khalatnikov perturbation theory based on the General 
Theory of Relativity combined with thermodynamics and a realistic equation of state for the pressure p = p{n,e), 
predicts that in the era after decoupling of matter and radiation small-scale inhomogeneities may grow very fast. As 
in Section IV, a fiat (i?(o) = 0) flrw universe with vanishing cosmological constant (A = 0) is considered. 



A. Introducing Observable Quantities 

The parameter /im (86) will be expressed in observable quantities, namely the present values of the background 
radiation temperature, r(o)^(tp), the Hubble parameter, Hitp) = cH{tp), and the redshift at decoupling, z(idoc)- 
The redshift z{t) as a function of the scale factor a{t) is given by 

m = - 1, (90) 

where a{tp) is the present value of the scale factor. For a flat flrw universe one may take a{tp) = 1. 
Substituting (78) into (86), one gets 



27r 1 /5 fcBr(o)(idoc) \ \ \ /ni\ 

Mm = T 7777 rVo 2 ' Adec = Aaidec, (91) 

where idee is the time when a perturbation starts to contract and Adcc the physical scale of a perturbation at time 
idee- Using (83) and (90), one finds 



2tt 1 J5 kBTfo-^^itp) 



= U f, A^il ' Adcc^Aa(tdec), (92) 
Adcc H(tp)|_2(tdoc) + IJ V 3 no- 
where it is used that T^g^{tdcc) = 2^(0)7 (^doc), and that T(o)^ oc after decoupling, as follows from (69). With (92) 
the parameter /im is expressed in observable quantities. 



B. Initial Values from WMAP 



The physical quantities measured by the Wilkinson Microwave Anisotropy Probe (wmap) [24] and needed in the 
parameter (92) of the generalised Lifshitz-Khalatnikov theory is the redshift at decoupling, the present values of 
the Hubble function and the background radiation temperature, the age of the universe and the fluctuations in the 
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background radiation temperature. The numerical values of these quantities are 

^(tdcc) = 1091, (93a) 

cH{tp) ^ n{tp) = 71.0km/sec/Mpc = 2.30 x 10"i^sec-\ (93b) 

T,o,^(ip) = 2.725 K, (93c) 

ip = 13.75 Gyr, (93d) 

<5T,(tdec) < 10-'^ (93e) 
Substituting the observed values (93a)-(93c) into (92), one finds 

15-69 

Mm = -r , Adoc m pc, (94) 

Adcc 

where it is used that m = toh = 1.6726 x lO^^^kg and 1 pc = 3.0857 x 10^^ m = 3.2616 ly. Boltzmann's constant is 
fee = 1.3806 X 10-23 J K-i. 

The WMAP observations of the fluctuations St (^dcc), (93e), in the background radiation temperature yield for the 
fluctuations in the energy density 

l'5.(idcc,g)| <10-^ (95) 

In addition, it is assumed that 

K{tdec,q)~0, (96) 

i.e., during the transition from the radiation-dominated era to the era after decoupling, perturbations in the energy 
density are approximately constant with respect to time. It follows from (81) that during the linear phase of the 
evolution one has Sn(t, q) « 5^{t^ q), so that the initial values 6n{tdcc, q) and (5n(tdcc, q) need not be specified. 



C. Star Formation in the Early Universe 



At the moment of decoupling of matter and radiation, photons could not ionise matter any more and the two 
constituents fell out of thermal equilibrium. As a consequence, the pressure drops from a very high radiation pressure 
p = IflBT^ just before decoupling to a very low gas pressure p = nk^T after decoupling. This fast and chaotic 
transition from a high pressure epoch to a very low pressure era may result in large relative pressure perturbations 
5p = pf]')/p(o). With (81) and (95) it follows that 5t = T^i)/Tio) could be large. As will be shown, relative initial 
pressure perturbations 

Spitdec, q) ~ Sritdcc, q) < -0.005, (97) 

may result in primordial stars, the so-called (hypothetical) population ill stars, several hundred million years after 
the Big Bang. 

The evolution equation (85) is solved numerically and the results are summarised in Figure 1, which is constructed 
as follows. For each choice of Sritdccq) equation (85) is integrated for a large number of values for the initial 
perturbation scale Adoc using the initial values (95) and (96). The integration starts at r = t/^doc = 1; i-G-j 
z{tdec) = 1091 and will be halted if either z = (i.e., r = [z(tdcc) + 1]'^^^), or 6^(1, q) — 1 for z > has been reached. 
One integration run yields one point on the curve for a particular choice of the scale Adcc if Se{t,q) = 1 has been 
reached for z > 0. If the integration halts at z{tp) = and still (5e(tp, q) < 1, then the perturbation belonging to that 
particular scale Adcc has not yet reached its non-linear phase today, i.e., at tp — 13.75 Gyr. On the other hand, if the 
integration is stopped at 5e(t, q) = 1 and z > 0, then the perturbation has become non-linear within 13.75 Gyr. This 
procedure has been performed for ST{tdcc,q) in the range —0.005, —0.01, —0.02, . . . , —0.1. Each curve denotes the 
time and scale for which 5e{t,q) = 1 for a particular 6T{tdcc, q)- 

The growth of a perturbation is governed by both gravity as well as heat loss. The smaller the scale of a density 
perturbation, the more it benefits from heat loss and the less from its internal gravity. From Figure 1 one may infer 
that the optimal scale for growth is around 6.2 pc ~ 20 ly. At this scale, which is independent of the initial value 
of the matter temperature perturbation STitdcci q), heat loss and gravity work together perfectly, resulting in a fast 
growth. Perturbations with scales smaller than 6.2 pc reach their non-linear phase at a later time, because their 
internal gravity is weaker than for large-scale perturbations. On the other hand, perturbations with scales larger than 
6.2 pc cool down slower because of their large scales, resulting also in a smaller growth rate. Since the growth rate 
decreases rapidly for perturbations with scales below 6.2 pc, this scale will be considered as the relativistic counterpart 
of the classical Jeans scale. The relativistic Jeans scale is much smaller than the horizon size at decoupling, given by 
dnitdcc) = 3rfdcc ~ 3.5 X 10^ pc « 1.1 x 10^ ly. 
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Star Formation starting at z = 1091 




10 20 30 40 50 

Perturbation Scale (parsec) at Decoupling 



Figure 1. The curves give the redshift and time when a linear perturbation in the energy density with initial values ^^(tdcc, q) ^ 
10~^ and 5e(tdcc,<7) ~ starting to grow at an initial redshift of z(tdcc) ~ 1091 becomes non-linear, i.e., 5^{t,q) — 1. The 
numbers at each of the curves are the initial relative perturbations SritdecQ) in the matter temperature. For each curve, the 
Jeans scale (i.e., the peak value) is at 6.2 pc. 



D. Heat Loss of a Density Perturbation during its Contraction 

The heat loss of a density perturbation during its contraction after decoupling can be calculated from the combined 
First and Second Law of Thermodynamics (56) rewritten in the form 

7^(0)4 = - — ('5n-^e)-^<5„, (98) 

where it is used that w = p^g^/e^o)- Substituting expressions (77) and (81) into (98) and using also (82), one finds the 
entropy per particle of a density perturbation: 

sfl{t,x)^ ^kB[3ST{to,x)^2Sn{t,x)], (99) 

where it is used that mc^ ^ fceTjo). From (99) it follows that for St < and ^„ > the entropy perturbation 
is negative, sfl^ < 0. Since for growing perturbations one has (5„ > the entropy perturbation decreases, i.e., 

sfl^ ^ — fcef^n < 0, during contraction. This is to be expected, since a local density perturbation is not isolated from 
its environment. Only for an isolated system the entropy never decreases. 



E. Relativistic Jeans Mass 



The Jeans mass at decoupling, -Mj(tdcc), can be estimated by assuming that a density perturbation has a spherical 
symmetry with diameter the relativistic Jeans scale Aj(idcc)- The relativistic Jeans mass at decoupling is then given 

by 

Mj(tdec) = y (5Adcc)^"-(o)(idcc)m. (100) 

The particle number density n(o)(tdcc) can be calculated from its value n^o^{tcq) at the end of the radiation-dominated 
era. By definition, at the end of the radiation-domination era the matter energy density n^„^mc^ equals the energy 
density of the radiation: 



(101) 
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Since n(o) oa a ^ and Tjo)-,. oc a one finds, using (90) and (101), for the Jeans mass (100) at time tdcc 

Mj(tdcc) = ^^AL"^%^ [^(^cq) + 1] [^(idcc) + 1]'. (102) 

Using (93a), the black body constant ae = 7.5657 x 10^^^ J/m'^/K'*, the red-shift at matter-radiation equaUty, 
z{toq) = 3196, WMAP [24], and the speed of hght c = 2.9979 x 10^ nr/s, one finds for the Jeans mass at decouphng 

Mj(tdoc)«3.5x lO^Mo, (103) 

where it is used that one solar mass IMq = 1.9889 x IQ-^" kg and the relativistic Jeans scale Aj(tdcc) = 6.2 pc, the 
peak value in Figure 1. 



VI. WHY THE STANDARD EVOLUTION EQUATION IS NO LONGER ADEQUATE 



The standard evolution equation for relative density perturbations S{t, x) in a flat, i?;,,) — 0, flrw universe with 
vanishing cosmological constant, A = 0, reads 



2H6 



■,2 



6^0. 



(104) 



As will become clear in the sequel, this equation holds true if and only if the pressure P(o) a-nd the energy density 
and their perturbations p^ij and £(1) obey the linear barotropic equation of state 



p — we. 



(105) 



The evolution equation (104) is derived from the Newtonian Theory of Gravity. In this derivation it is assumed that 



ivs/cy 



w 



in the radiation-dominated era, and w — and Vs/c is time dependent and given by (78) after 



decoupling of matter and radiation. 

The fact that the factor of proportionality w is a constant can be demonstrated using the general expressions (29) 
and (44). From (11) and (105) one finds that ^ w and Pn = 0. This implies with (29) that f3'^ = Pe = w, so that 
with (44), one gets w = 0. Thus, thermodynamics and the conservation laws (7c) and (7e) imply that for an equation 
of state p = we, the factor of proportionality w must be constant. 

In order to compare equation (104) with the generalised Lifshitz-Khalatnikov theory (43), the relativistic analogue of 
(104) pertaining to the equation of state (105) will now be derived from the background equations (7) and perturbation 
equations (40), with — and A = 0. Since, by (105), one has p„ = 0, equations (7e) and (40d) need not be 
considered. With S = e^i^/e^o), equation (40b) becomes 



S+{l + w) 



-^(1) 
4H 



0, 



(106) 



where kEjoj = 3H^, (7a), has been used. Differentiating (106) with respect to time and eliminating the time-derivatives 
of H, and with the help of the background equations (7) and perturbation equations (40), respectively, and, 
subsequently, eliminating with the help of (106), one finds that the system of equations (40) reduces to 



S + 2H5 - 



;Ke 



+ i/(2 - 3u;)i5(i, + 



iw) 



S = -3Hw{l + w)l9^,y, 



w 

1 + w a? 



= 0, 



(107a) 
(107b) 



where it is used that w — Q. The system (107) consists of two relativistic equations for two unknown quantities, 
namely the density fiuctuation 5 and the divergence of the spatial part of the fluid four-velocity. Thus, the 
relativistic perturbation equations (40) which are derived for a general equation of state for the pressure p = p(n, e) 
reduce to the relativistic system (107) for the equation of state (105) with constant w. 

In the derivation of equations (40) no assumptions or approximations (other than linearisation of the Einstein 
equations and conservation laws) have been made. Since the relativistic system (107) is derived from the system 
(40), the set (107) is exact and is, therefore, superior to the standard Newtonian equation (104) which is only an 
approximation. As a consequence, equation (104) is only adapted to the equation of state (105) and the coefficient 
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(vs/c)^ must be replaced by the constant w. In other words, the equation of state p = we requires that the speed of 
sound Vs{t) is for all times equal to its initial value Vs{to) = c^/w. 

The gauge modes (38) are solutions of the set (40), implying that the gauge modes 

^(t,a.) = '^MmW = -3i7(i)VX^)(l + ^), ^<.,(t,^) = (108) 

are solutions of equations (107), whereas (5 is a solution of (104) for constant V', see (64). Therefore, equations (104) 
and (107) have no physical significance, since the unknown gauge function '^(x) cannot be fixed by imposing physical 
initial conditions, as has been demonstrated in Section IIIH. 

Finally, the manifestly covariant and gauge-invariant analogue of equations (107) will be given. For i?(o) = 0, A = 
and an equation of state p — we with p„ = and constant {3"^ = w, the set (43) reduces to 



(5e + H{1 - Qw)6e - 



+ ^K£(o)(l - 9w)(l + 3w) 



5e = 0, (109a) 



1 d_ 



On - 



1 



= 0. (109b) 



For w = |, these equations become equal to (68). However, for w = Q the set (109) differs substantially from the 
set (79). This is due to the fact that, in (109), (P = w = Q and p = Q. In contrast, in equations (79) p 7^ and 
P{t) « Vs{t)/c, (78). As a consequence, the term $/f3 = —H is included in (79), whereas it is missing in equations 
(109). 

The sets (107) and (109) are both derived from the basic perturbation equations (40). The fact that the sets (107) 
and (109) differ considerably from each other is that (107) is based on the gauge-dependent variable 6 = e(i)/e(o), 
whereas (109) is the evolution equation for the gauge-invariant quantity Sg = ef',/e(o), so that 5^ ^ 6. As follows from 
(37) with S = e and (64), £(i) is gauge-dependent also in the non-relativistic limit. That is why efl^ will never become 
equal to e^iy 

In both sets (107) and (109) the divergence of the spatial part of the fluid four-velocity is taken into account, 
whereas is missing in the standard equation (104). 

The differences and similarities between equations (104), (107) and (109) will be discussed in the next two subsec- 
tions. The conclusion will be that the standard equation (104) is inadequate throughout the history of the universe, 
whereas (107) and (109) are inadequate only after decoupling. 



A. Radiation-dominated Era 



In this era the fluid consists of a mixture of radiation and matter. Since radiation gives by far the largest contribution 
to the pressure, the equation of state for the pressure is to a very good approximation given by p = ^e (i.e., fP = w = ^) 
so that pn = 0, implying that the right-hand side of equation (43a) vanishes and that this equation reduces to (109a). 
Since from the outset a realistic equation of state for the pressure p = p(n, e) has been used in the derivation of the 
generalised Lifshitz-Khalatnikov theory, an additional equation (43b) is found. This equation survives in the form 
(109b). As a consequence, the generalised Lifshitz-Khalatnikov theory predicts that perturbations in the particle 
number density are gravitationally coupled to perturbations in the radiation density, irrespective of the nature of the 
particles. This may prevent CDM from contracting faster than ordinary matter. This phenomenon can neither be 
discovered by the standard perturbation equation (104) nor by the exact equations (107), because these equations are 
adapted only to the equation of state (105), without taking into account the particle number density. Therefore, the 
equation of state (105) is incomplete. 

For large-scale perturbations V^J — >■ both (104) and the homogeneous part of (107a) have as solutions S oc t and 
S cx t~^, where the latter is the gauge mode (108). For V^i5 — )■ the solution of equation (107b) is the physical mode 
cx cx t~^/^, implying that the corresponding particular solution of (107a) is the physical mode S cx i^/^. The 
physical modes of the system (107), i.e., S (xt and S cx i^^^, are precisely the two parts of the general solution (75) of 
the manifestly covariant and gauge-invariant equation (109a). This shows that the physical part of is contained 
in equation (109a). Since -d^i-, contributes only to the slow growing mode 5^ cx t^^^, plays only a minor role in 
large-scale perturbations. That is why the new perturbation theory (109) with i^jij 7^ predicts the same growth rate 
5e oc t as the standard theory (104) which does not take into account. 

If is put to zero in equation (107a) one arrives at the standard equation (104). This is not allowed, since 
■iJji) is an essential quantity for the evolution of density perturbations. Indeed, if there is no fluid flow towards a 
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density perturbation, then that perturbation cannot grow. The generahsed Lifshitz-Khalatnikov perturbation theory 
(109), which has ^ 0, predicts that small-scale density perturbations, for which V^(5 is large, oscillate with an 
increasing (76) amplitude. In contrast, the standard equation (104), which has = 0, predicts that small-scale 
density perturbations oscillate with a constant amplitude. Thus, as can be inferred from (107), ^^(i) plays a major 
role in the evolution of small-scale density perturbations. It is, therefore, important that the divergence of the spatial 
part of the fluid four- velocity is taken into account. Since = in the standard equation (104), this equation 
is inadequate to study density perturbations in the radiation-dominated era of the universe. 



B. Era after Decoupling of Matter and Radiation 

In this era one has w <C 1, so that (104) and (107) yield nearly the same results, since t9(i) plays only a minor 
role. In fact, in the limit w — >■ 0, the solution of (107b), -iJji) cx a~^, is equal to the gauge mode (108), implying 
that the physical part of physical, is exactly zero. This is to be expected since p — ?> implies physical — >■ 0, 

and vice versa, see the text below (61). Therefore, equation (107b) need not be taken into account and "iJji), being 
only a gauge mode, may be put to zero. In this case, the solutions of (104) and (107a) are i5 oc t"^^^ and 6 oc t~^, 
where the latter is the gauge mode (108) with constant ?/'• Since the sets (107) and (109) are both derived from the 
basic perturbation equations (40), the physical part of "iJji) vanishes also in the set (109) for w — > 0. Therefore, the 
sets (107) and (109) should yield, apart from the gauge modes, exactly the same solution. Indeed, the solutions of 
(109a) are the modes 5^ oc t"^^^ and 5g cx Ignoring the decaying solutions, the three equations (104), (107) and 

(109) yield exactly the same result, namely growing density perturbations with a scale-independent growth rate. This 
result cannot be correct for the following two reasons. Firstly, since there is no fluid flow, i.e., physical — >■ 0, in the 
limit p ^ 0, perturbations do not grow and the gravitational field is static (65b). Secondly, given two energy density 
perturbations with the same energy density £(0) -|- e^',, then the large-scale perturbation has a stronger gravitational 
field than a small-scale perturbation. As a consequence, growth rates are scale-dependent. The incorrect results of 
equations (104), (107) and (109) are due to the equation of state (105), which is incomplete and, therefore, does not 
describe the cosmic fluid adequately in the perturbed universe after decoupling. 

Since the generalised Lifshitz-Khalatnikov theory (43) is adapted to the realistic, non-barotropic, equation of state 
for the pressure p = p{n, e), it does not have the above mentioned disadvantages of the standard perturbation theory 
(104) and equations (107) and (109). Firstly, the equations of state (77) yield for the quantity /?, (29), the correct, 
time-dependent, speed of sound (78), whereas = w is constant for an equation of state p = we. Secondly, it follows 
from Figure 1 that in the new perturbation theory the growth rate depends on the scale of a perturbation, as it should 
be. There is a sharp lower limit, namely 6.2 pc, below which density perturbations do not grow at all (89). 
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##################################################################################### 

> # Relativistic Cosmological Perturbation Theory and the Evolution of Small-Scale Perturbations 
^ # 

_> # P . G . Miedema 

_> ##################################################################################### 
_> # Maple 16 

=^ # 

_> # Consistency check of the perturbation equations (43) . 
=?■ * 

_> # To check Eqs.(43) with coefficients (45), substitute (42) into Eqs.(43) 
_> # and use the background equations (7) and the perturbation equations (40) 
_> # to evaluate all time derivatives . 

_> ##################################################################################### 

_> restart; 

_> # Constants : 

_> # kappa = 8*pi*G/c'^4, 

_> # Delta is the nabla operator, 

I> ##################################################################################### 
_> # Coefficients b_l, b_2, b_3, (45) : 

> bl : =kappa*e_0 (t) * (l+w(t) ) /H (t) -2*dif f (beta (t) ,t) /beta (t ) -H (t) * (2+6*w (t) +3*beta (t) ^2 ) +R_0 (t)* (1/(3*H 
(t) )+2*H(t) *(l+3*beta(t) ^2) / (R_0 (t) +3*kappa*e_0 (t)* (l+w(t) ) ) ) ; 



bl : 



Ke_OU) {1 +w(t)) 



P(0 



-W) {2+6w{t)+3^[t) ) +R_0{t) 



1 



2 Hit) (1 +3|3W^) 



Hit) p(0 ''[3H(.t) R_0{t)+3Ke_0(t)(l+w{t)) 

~> b2:=-(l/2) *kappa*e„0 (t) * (l+w(t) ) * (l+3*w(t) )+H(t) ^2* (l-3*w (t) +6*beta (t) ^2* (2+3*w(t) ) )+6*H(t) * (diff 
(beta (t) , t) /beta (t) ) * (w (t) +kappa*e_0 (t) * (1+w (t) ) / (R_0 (t) +3*kappa*e_0 (t) * (1+w (t) ) ) ) -R_0 (t) * ( (1/2) *w 
(t) +H (t) ^2* (l+6*w (t) ) * {l+3*beta (t) ^2) / (R_0 (t) +3*kappa*e_0 (t) * (1+w (t) ) ) ) -beta (t) ''2* (Delta/a (t) ^2- 
(l/2)*R_0(t)); 



bl--- - Ke_0{t) (1 +w{t)) {1+3 w{t)) +mf (1 -3 wit) +6 P(fl (2+3 w{t))) 



6m)(^m)] Mt) + — ^ 

{dt J [ R 0(f) +3 K 



Ke Oit) (1 +w(0) 



e Oit) (1 +wit)) 



P(f) 



-« Oit) I 1 Mt) + //(0^(l+6w(f))(l+3P(f)^) 
2 R_0{t) +3 K.e_Oit) (1 +w{t)) 



- m 



ait)^ 



RJKt) 



> b3 : = (-18*H (t) ^2* (e_0 (t) *p_en (t) * (1+w (t) ) +2*p_n (t) *dif f (beta (t) , t) /beta (t) / (3*H (t) ) +p_n (t) * (p_e (t) ■ 
beta (t) ^2) +n_0 (t) *p_nn (t) ) / (R_0 (t) +3*kappa*e_0 (t) * (1+w (t) ) ) +p_n (t) ) * (Delta/a (t) ^2- (1/2) *R_0 (t) ) * 
(n_0(t)/e_0 (t)); 



b3 : 



ejit) 



18 //(fl 



ejit) p_enit) (1 +w(/)) + 



2 

3 p(fl Hit) 



+p_nit) [p_eit) - P(O^) +n_Oit) p nnit) 



R_0it)+3 Ke_Oit) (1 +wit)) 



+p_n{t) 



ait)-" 



ROit) n_Oit) 



J 

_> ##################################################################################### 
_> # Background equations 

_> ##################################################################################### 



> w(t) :=p_0(t)/e_0(t); 



_> # Scale factor 

_> define (a, diff (a (t) , t) =H (t) *a (t) ) , 
> diff (a(t) ,t) ; 



w(f) : 



. P_Oit) 
e_Oit) 



Hit) ait) 



> # The time derivative of Eq. (7a) combined with Eqs . (7b) and (7c) yields the 
_> # time derivative of the Hubble function H: 

~> define (H, diff (H(t) ,t)=- (1/6) *R_0 (t) - (1/2) *kappa* (e_0(t)+p_0 (t) ) ) ; 

> diff (H(t) ,t) ; 



- - R_0{t) - - K (e_0(0 +p_0{t)) (6) 
_ 6 2 

_> # Either define the three-space curvature R_0 by Eq. (7b) : 

_> define (R_0, dif f (R„0 (t) ,t) =-2*H(t) *R_0 (t) ) ; 

_> # or by Eq. (8) with k=+l, 0, -1: 

_> # k:=+l; R_0 (t) :=-6*k/a(t) ^2; 

~> diff (R_0(t) ,t) ; 

_ -2H(t)R_0{t) (7) 

_> # Energy conservation law (7c) : 

_> define (e_0, diff (e_0(t) ,t)=-3*H(t) *e_0 (t)*(l+w(t) ) ) ; 

~> diff (e_0(t) ,t) ; 

# Particle number conservation law (7e) : 

_> define (n_0, diff (n_0(t) , t) =-3*H (t) *n_0 (t) ) ; 

~> diff (n_0(t) ,t) ; 

_ -3H{t)n_0U) (9) 

_> # Time derivative of p„0 (t) : 

_> define (p_0, dif f (p„0 (t) ,t) =p_e(t) *diff (e_0 (t) ,t) +p_n(t) *diff (n_0 (t) ,t) ) ; 

> diff (p_0(t) ,t) ; 

-3 p_e{t) Hit) ej){t) [ 1 + EJI^ 1 - 3 p_n(t) H(t) nj)[t) (1 0) 

- I ej)(t) J 

_> # Partial derivative p_e(t) of the pressure p: 

_> define (p_e, diff (p_e(t) , t) =p_ee (t) *dif f (e_0 (t) , t) +p_en (t) *dif f (n_0 (t) ,t)); 

> diff (p_e(t) ,t) ; 

-3 p_ee(t) H(t) e_0(t) [ 1 + EJI^ \ _ 3 p_en{t) H{t) n_0(t) (11) 
_ I e_0{t) J 

_> # Partial derivative p_n(t) of the pressure p: 

_> define (p_n, diff (p_n(t) , t) =p_en (t) *dif f (e_0 (t) , t) +p_nn (t) *dif f (n_0 (t) ,t)); 

~> diff (p„n(t) ,t) ; 



-3p en(t)H[t)e 0(0 f 1 + 1 - 3 n nn{t) H{t) n 0(t) (12) 

- y e 0(t) J ^ ' 



# Quantity beta(t) : 

beta(t) :=sqrt (diff (p_0 (t) ,t)/diff (e_0 (t) ,t) ) ; 



P(0 : 



-9 p e(t) HU) e 0(t) [ 1 + EJI}^ \-9p n{t) Hit) n 0{t) 
- I e 0(t) ) - 

p_0(t) 



Hit) e_0(t) 1 + 



e_0(t) 



(13) 



_> ##################################################################################### 

^ # First order perturbation equations 

_> ##################################################################################### 
_> # First order perturbation to the pressure: 

> p_l(t) :=p_e(t)*e_l(t)+p_n(t)*n_l(t) ; .... 
_ pj(t)--p_e(t)ej(t)+p_n(t)nj(t) \\'\) 
J> # Energy conservation law (40b) : 

> define (e_l, diff (e_l (t) ,t)=-3*H{t) * (e^l (t)+p_l (t) )-e_0(t) * (l+w{t) ) * (theta (t) + (kappa*e„l (t) + (1/2) *R_1 
_ (t))/(2*H(t)))); 

> diff (e_l (t) ,t) ; 



-3 H(t) (ej(t) +p_e{t) ej(t) +p_n(t) nj{t)) - e_0(t) 1 + 



pjm_ 

e_0[t) 



m + 



, Ke l(t) + - R l(t) 

1 2 ~ 

2 H(t) 



(15) 



J> # Particle number conservation law (40d) : 

_> define (n J, diff (n_l (t) ,t)=-3*H(t) *n_l (t)-n_0 (t) * (theta (t) + (kappa*e_l (t) + (l/2) *R_1 (t) ) / (2*H(t) ) ) ) ; 
~> diff (n_l (t) ,t) ; 



-3 H(t) n_l(t) -n_0(t) 



9(0 + 



, Ke ht] + - R Ht) 
1 - 2 - 



2 Hit) 
_> # Momentum conservation law (40c) : 
> define (theta, dif f (theta (t) , t)=-H(t) * (2-3*beta(t) ^2) *theta (t) -Delta/a (t) ^2*p„l (t) / (e„0 (t) * (l+w(t) ) ) ) 



(16) 



> dif f (theta (t) ,t) ; 



-Hit) 



2 + 



, -9p e(t) Hit) e Oit) (\+ EJ!MI \-9p nit) H[t) n 0(f) 
1 ^- - \ e Oit) ' 



Hit) e_Oit) 1 + 



pjm_ 

ejit) 



0(0- 



A ip_eit) eJit) +p_nit) njit)) 



ait)^e_Oit) 1 + 



pjm_ 

e_Oit) 



(17) 



_> # Evolution equation for the local perturbation to the 
_> # global spatial curvature, (40a) : 

> define (R_l, dif f(R_l (t) ,t)=-2* (H(t) *R_1 (t) -kappa*e_0 (t) * (1+w (t) ) *theta (t) ) -R_0 (t) / (3*H(t) ) * (kappa* 
_ e_l(t) + (l/2)*R_l(t))) ; 

> diff (R_l (t) ,t) ; 

„ 1 R^Oit) Uejit) + J R_m 

-2Hit)R lit) +2Ke Oit) 11 + | 6(0 

_ - - I e_Oit) ) 3 Hit) 

_> ##################################################################################### 

_> # Gauge-invariant density perturbations 

_> ##################################################################################### 
_> # Gauge-invariant perturbation to the energy density, (41a) : 

> e_gi (t) : = (e_l (t) *R_0 (t) -3*e_0 (t) * (1+w (t) ) * (2*H (t) *theta (t) + (1/2) *R_1 (t) ) ) / {R_0 (t) +3*k;appa*e_0 (t) * 
(l+w(t))); 



(18) 



e_giit) : 



e lit) R Oit) -3e Oit) I 1 + ] ( 2 Hit) 6(0 + - R Ht) 

- - - ' e Oit) J [ 2 - 



R Oit) +3 Ke_0(O 1 + 



pjit) 
e_Oit) 



(19) 



_> # Gauge-invariant perturbation to the particle number density, (41b) : 
> n„gi (t) : =n_l (t) -3*n_0 (t) * (kappa*e_l (t) +2*H (t) *theta (t) + (1/2) *R_1 (t) ) / (R_0 (t) +3*kappa*e_0 (t) * (1+w (t) 

)); 



n_giit) :=«__/ (0 ■ 



3 n_0(0 I Kejit) +2 Hit) 6(0 + - RJit) 



R_0it)+3 Ke_Oit) 1 + 



pjm_ 

e_Oit) 



(20) 



_> # Gauge-invariant contrast functions, (42) : 
> delta_e (t) : =e_gi (t) /e_0 (t) ; delta_n (t) : =n_gi (t) /n_0 (t) ; 



delta eit) := 



e !it)R 0(0 -3 e 0(0 fl +^=^1 f 2 //(O 6(0 + - « Aol 

~ ~ ~ I e_0(0 J \ 2 - ) 

R Oit) +3 Ke_0(0 f 1 + f^J^ 1 \ eJit) 



nJit) ■ 



e_Oit) 

3 nJit) I Ke_/(0 +2 Hit) 6(0 + ^ RJit) 



R Oit) +3 Ke 0(0 1 + 



delta_nit) := ■ 



pjm_ 

e Pit) 



n_Oit) 



(21) 



_> # Consistency check for Eq. (43b) [left-hand side minus right-hand side should yield zero] : 
> simplify (diff (delta_n(t)-delta_e(t) / (l+w(t) ) ,t)-3*H(t) *n_0 (t)*p_n(t) / (e_0{t) * (l+w(t) ) ) * (delta_n(t)- 
delta„e(t)/ (l+w(t) ) ) ) ; 



(22) 

_> # Consistency check for Eq. (43a) [left-hand side minus right-hand side should yield zero] : 
> simplify (diff (delta_e (t) , t$2) +bl*dif f (delta_e (t) , t) +b2*delta_e (t) -b3* (delta_n (t) -delta_e (t) / (1+w (t) 



))); 

# Conclusion: the set (43) is correct. 



(23) 



27 



